Signatures of surface states in bismuth at high magnetic fields 
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Electrons in a metal subject to magnetic field commonly exhibit oscillatory behavior as the field 
strength varies, with a period set by the area of quantized electronic orbits. Recent experiments 
on elemental bismuth have revealed oscillations for fields above 9 tesla that do not follow this 
simple dependence and have been interpreted as a signature of electron fractionalization in the 
bulk. We argue instead that a simple explanation in terms of the surface states of bismuth exists 
when additional features of the experiment are included. These surface electrons are known to have 
significant spin-orbit interaction. We show the observed oscillations are in quantitative agreement 
with the surface theory, which we propose to test by studying the effect of the Zeeman coupling in 
higher fields, dependence on the field orientation, and the thickness of the samples. 



Elemental bismuth has recently attracted renewed in- 
terest due to the experimental observation by Behnia et 
al. of anomalous quantum oscillations at high magnetic 
fields > 9 T along the trigonal crystal axis |l| . Plateau- 
like features in the Hall resistivity in this range of fields 
were taken as an indication of bulk electron fractionaliza- 
tion, in a manner reminiscent of the fractional quantum 
Hall (FQH) effect in two dimensions (2d). A second ex- 
periment by Li et al. found evidence for a new phase of 
electrons at such high fields [2j. Some theoretical work 
has followed 0, however, the anomalies [l[ remain 
unresolved. 

The idea of fractionalization in a bulk material is quite 
intriguing, especially in an isotropic semimetal such as 
bismuth. To this date, the only known realizations of 
fractionally charged particles are solitonic excitations of 
conducting poiymersjfjj] , FQH state of 2d semiconductor 
heterostructures 0, 0], and confined quarks of quantum 
electrodynamics. Highly anisotropic relatives of the inte- 
ger quantum Hall effect are found in layered semiconduc- 
tor heterostrctures § and Bechgaard salts 0, [l(| but no 
FQH state yet. The key experimental signature of FQH 
states is a plateau in the sideways resistivity at fractional 
values of h/e 2 and a concomitant vanishing of the longi- 
tudinal resistivity. While the stability of FQH states de- 
rives from Coulomb interaction and weak disorder, their 
existence owes to the essential role of topology in 2d. 
Theoretical extensions of fractionalization in 3d rely ei- 
ther on an anisotropic structure or on background topo- 
logical structures such as solitons in a relativistic Dirac 
Hamiltonian. It is not clear that such anisotropy as in the 
first scenario occurs in bismuth and although electrons in 
bismuth are known to have a Dirac-like dispersion, it is 
not clear what field would provide the topological back- 
ground in the second one. Furthermore, the longitudinal 
resistivity in bismuth is found to be rather featureless at 
such field orientations [H, [ill ] . The torque magnetome- 
try which is a sensitive bulk measurement, does not 
show these anomalies either, casting further doubt on a 
FQH scenario. We believe based on these considerations 
an alternative explanation of the observed anomalies is 
strongly favored. 



Indeed we shall show that they can be explained in a 
simple fashion by the states confined to the (111) surface 
of bismuth. In particular, a distinct indexing of Landau 
levels according to the surface states with a period of 
oscillation w 0.016 T -1 explains the experimental data. 
The surface theory predicts: (1) the existence of addi- 
tional features that might have already been observed 
and an additional peak at higher fields ~ 60 T corre- 
sponding to the surface quantum limit; (2) a distinct 
dependence of the field at the peaks on the angle 9 of a 
tilting field relative to the trigonal axis. The magnetic 
fields so far studied fall below the surface quantum limit. 
However, the Zeeman coupling of surface electrons could 
reduce the quantum limit down to 40-50 T. These fea- 
tures should allow for a falsification of the surface-state 
theory of the anomalous peaks. 

Quantum oscillations are ubiquitous in metals. As the 
field changes, the Landau levels cross the Fermi energy 
and depopulate their electrons at certain fields B~ 1 = 
(n + 7)A£?~ 1 , where < 7 < 1 is a correction arising 
from the quantum mechanical nature of the electronic 
orbit (7 = 1/2 for free electrons) and 



Here, Sp is an extremal area of the Fermi surface per- 
pendicular to the field and $0 = hc/e is the unit flux 
quantum. The coincidence of the Fermi and Landau en- 
ergies results in enhanced or singular contributions to 
most electronic properties; hence they show an oscilla- 
tory structure with B^ 1 with the period of oscillations 
given by the Onsager relation, Eq. fl}. 

Due to its low carrier density and long mean-free path, 
bismuth exhibits periodic oscillations in its susceptibility 
and resistivity starting at fields below 1 T [T^, [HJ • Us- 
ing Eq. (fl]) and the known Fermi surfaces of electron- 
and hole-like carriers in bismuth, one finds that 9 T cor- 
responds to depopulating all but the lowest Landau level 
for holes, followed closely by that of electrons, the so- 
called quantum limit. A recent study of the Landau levels 
of electrons and holes in bismuth, including the effects of 
the linear dispersion of electrons, accounts for most of the 
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peaks observed below this field Q with distinct depen- 
dence on the field orientation. The anomalous peaks [lj 
do not follow Eq. (JXJ) . 

Bismuth also has a variety of interesting electronic 



states confined to its surfaces [14|. The (111) surface per- 
pendicular to the trigonal axis has been studied in depth 
using ARPES and other techniques and found to host a 
number of carriers equivalent to a 1 /im-thick sample of 
the bulk. The samples in the experiment [l| have a thick- 
ness of 0.8 mm; so, despite their higher relative carrier 
density one expects the Nernst signal of the surface to be 
weaker than that of the bulk by a factor of at least 1000. 
This is certainly true for fields below the quantum limit. 
However, above the quantum limit, the strong peaks from 
the bulk will be absent and the surface signal could be 
detected much more easily. The surface states in bis- 
muth and related materials have also gained attention in 
relation to "topological insulators" [15j. The important 
ingredient in such physics is the significant spin-orbit in- 
teraction (SOI) that causes a large energy splitting in 
the surface states and leads to a single spin state per 
momentum. In bismuth the role of SOI has long been 
appreciated and has been shown to_be the cause for the 

At the (111) sur- 
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existence of hole Fermi surfaces 
face, the SOI results in six elongated hole Fermi pockets 
arranged around an electron pocket centered at the T 
point in the hexagonal Brillouin zone. 

In order to study the effects of the Zeeman coupling 
and tilting fields we model the surface states of Bi(lll) 
by a low-energy continuum hamiltonian 
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where p is the momentum operator, m is the effective 
mass of electrons, <r = (a x ,a y ) are spin Pauli matrices, 
and the last term is the Rashba-type SOI with strength a. 
The trigonal axis is taken to be along the z direction. The 
spectrum of Eq. @ is shown in Fig. [TJA-B and consists of 
a small and a large Fermi pocket. There is a single spin 
state for each momentum p given by a spin along pxz. In 
this model both of these pockets are electron-like, which 
is an artifact of the model. In reality, for higher momenta 
away from the T point, the lattice effects become impor- 
tant. This results in a bending of the bands, giving rise to 
hole-like pockets. The lattice anisotropy further reduces 
the symmetry creating disjoint oblongular hole pock- 
ets. Based on angle-resolved photoemission (ARPES) 
data [3], we set m = 0.45m e (m e is the real electron 
mass), a = 1.2 eVA and tp = 68 meV. These param- 
eters reproduce the radius of the small electron pocket 
p e F = -nia+^f{ma) 2 + 2me F = 0.043 A" 1 (see Fig. HP) 
within the experimental resolution (~ 0.023 A -1 ) and 
give a Fermi velocity vf = Ppjrn + a = 1.9 eVA in good 
agreement with the data. The lower band crosses the 



Fermi energy at p F = p e F + 2ma = 0.185 A 
sponding to one end of the hole pocket. 
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FIG. 1: (color online) The spectrum of the surface hamilto- 
nian (A) without and (B) without spin-obit interaction. The 
Fermi energies are marked, as are the Fermi wavevectors for 
the small and large Fermi pockets. (C) Sketch of Fermi pock- 
ets in the surface Brillouin zone (to scale). The two circles 
are from our model. The central (yellow) hexagon and the 
outward (green) ellipses are from ARPES measurements [l4l ]. 
The thickness of the lines (pink) represents the experimental 
resolution. The dashed hexagon indicates the symmetry of 
the pockets and the Brillouin zone. 



In 2d, the period of oscillations is also given by Eq. 
but with Sf now simply being the area of the Fermi 
pocket. We shall see this explicitly in our model. The 
orbital effect of the magnetic field is found by the minimal 
coupling p — > p — | A in Eq. ^ , where A is the vector po- 
tential. The Zeeman energy is Hz — —gfiB&zB for a field 
along z, where fin — eh/(2m e c) is the Bohr magneton. 
The experimental determination of the g-factor usually 
relies on matching the Zeeman energy splitting with the 
Landau level spacing, and thus requires a knowledge of 
the spectrum itself. In the following we shall take g as 
a fitting parameter which allows its determination from 
our theory in a tilting field. 

The spectrum in the magnetic field can be found ex- 
actly. Let us first define 
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corre- with <7 c ff = g — m e /m. Then by using the the commuta- 
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We have assumed the spatial variation of the field to 
be negligible, which is justified when the last term in 
Eq. (T4| is small compared to the Landau level spacing 
~ hcu c . This yields SB/B <C L/X where SB is the spatial 
variation of the field over a characteristic length L, and 
A = .g c fffi/(4m e a) = g e s x 1.6 A. This is in fact a much 
less strict condition than the one needed for obtaining 
Landau levels when g c g = 0, namely, 5B/B n < n for 
the nth Landau level. 

The spectrum of II is the same as that of Dirac elec- 
trons in a magnetic field with velocity a and Zeeman 
coupling <? c ff . This is known in the context of the quan- 
tum Hall effect in graphene, and is given by 



E n ,s — ± 



J 2mhu c (n + s) + (^^B f. (6) 



Here n — 0, 1, • • • and s = 0, 1 are Landau level indices 
and lo c = eB / (mc) is the cyclotron frequency. From 
Eq. ([5]), we have 



e^ jS = hw c {n+s)±y/2ma 2 huj c {n + s) + (g c sfiBB) 2 . (7) 

The ± correspond to small (— ) and large (+) Fermi pock- 
ets. Note that the level n, s = 1 are degenerate with 
n + l,s = 0. The lowest Landau level n — s = is 
non-degenrate. 

Let us first neglect the effective Zeeman coupling by 
setting g c ff = 0. Then by setting e„ iS = tp we can derive 
Eq. {T]) with Sf — Sp = Tr(pp) 2 for the small electron 
pocket and similarly for the large one. Consequently, the 
fundmamental period of oscillations in the surface theory 
is then found to be 
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The Zeeman term changes the linear dependence of 1/ B n 
on the Landau level index. This can be understood 
by formally neglecting the orbital coupling and keep- 
ing only the Zeeman coupling Hz- Then one finds 
two bands e(p, a) — §^ + o^Ja 2 \> 2 + (g cS ^BB) 2 , where 
u = ±1 is the sign of the spin projected along ap x 
z + g ff/!Bi?z. As a result, the small electron pocket 
shrinks and the large Fermi pocket expands by the same 
area AS F = 4^ i^^-B) 2 . The Fermi energy is un- 

Pf+Pf a ' 

changed. With our choice of parameters, ASf/S f = 
7.4 x 10- 7 {g cff B/T) 2 w 10% when g cS B = 370 T. The 
field-dependent area of the Fermi pocket then leads to a 



nonuniform Landau level spacing and the depopulation 
of Landau levels at lower fields. In particular the sur- 
face quantum limit -Bi(g ff) is a decreasing function of 
<? c ff . The deviation from the linear dependence should be 
observable in fields > 40-50 T for a value of g e g ~ 10. 

We now turn to the experiments. Ref. [l| identified three 
clear peaks in the Nernst signal beyond the bulk quan- 
tum limit. In Fig. [2)3, we fit these peaks with Eq. 0. 
The index n now counts the Landau levels of the surface. 
The quantum limit is found to be B\{g c g — 0) = 63 T 
when the Zeeman coupling is neglected. Also shown are 
fits obtained for g e g = 8 and 15 with B\ = 54 T and 
43 T, respectively. In this surface scheme, the peaks in 
Ref. Ill correspond to n = 2, 3 and 5. The bulk quantum 
limit coincides with the surface n — 7 Landau level. The 
n = 4 and 6 Landau levels correspond to fields 15.7 T 
and 10.5 T, respectively. We reproduce in Fig. the 
data of Ref. [l| with the above features marked. Interest- 
ingly, the n — 4 and 6 Landau levels seem to show up in 
this data, in addition to the more clear peaks originally 
identified. The n — 6 level is strongly shadowed by the 
bulk quantum limit signal, and the n = 4 peak seems to 
develop for temperatures < 0.83 K and is clearly visible, 
though rather broad, at 0.56 K. 

We submit that the good fit obtained here with the 
extended set of peaks observed in Ref. [H is evidence of 
the surface origin of these peaks. More recently, further 
magnetotransport experiments by Huber et al. (l6| on 
bismuth nanowires in fields up to 14 T have also clearly 
identified peaks caused by surface Landau level structure. 
The period of oscillations in nanowires is ~ 0.025 T _1 , 
larger than that found here. This could be caused by the 
different geometry of the nanowire surfaces which can 
influence the size of the Fermi pockets. 

We briefly discuss the effects of tilting the magnetic 
field on the surface states. For a field at an angle 6 with 
the trigonal axis, the effective field entering the orbital 
coupling is reduced by a factor cos 9. The Zeeman cou- 
pling should also include a term — <7_L/iscr • B^, where 
Bj^ is the in-plane component of the field. That the 
orbital coupling is not affected by B^ can be seen by 
choosing a gauge where A = (0, Bx cos 8, |B_l x r|) where 
r = [x, y) is the in-plane coordinate vector. The in-plane 
Zeeman coupling has the effect of shifting the momenta 
along Bj^ x z. In addition, it results in an increase in 
the Fermi energy and introduces an angular dependence 
in the spectrum thus making the Fermi surface slightly 
anisotropic with a remaining reflection symmetry around 
Bj^. However, unlike the z-axis Zeeman coupling, it does 
not lead to a significant change in the size of the Fermi 
pockets. Its strength is further diminished at small 9 
and therefore we do not expect a significant effect from 
the in-plane field on the surface states and the period of 
quantum oscillations, unless g± is anomalously large. 

Shubnikov-de Haas measurements of single-crystal 
thin films (<~ 10 fim) of bismuth up to a few Tesla show 
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field is decreased relative to the bulk. 

There are additional electron pockets near the M point 
(not shown in Fig.[T]C) which could give rise to additional 
oscillations. The holes are nearly compensated with the 
electrons. Therefore one should expect a beating with 
a period < 6 times the period of the central electron 
pocket. The resulting nodes in the amplitude of oscilla- 
tions might offer an explanation for the suppression of 
the peak near n = 4. A clear test of the surface origin 
of the anomalous peaks is the dependence on the thick- 
ness of the samples. One could also try contacting only 
the surfaces for comparison. A whole family of bismuth- 
based materials show similar surface states offering an- 
other venue for testing their transport signatures. This 
is especially illuminating in topological insulators [IH 
where the bulk states are gapped. It is conceivable that 
the anomalies are caused by a bulk reorganization of elec- 
trons [§[ . This is even more plausible if accompanied by 
a drop in the longitudinal resistivity. However, absent 
further evidence, we believe the surface theory provides 
a satisfactory explanation of the current data. 

The authors acknowledge useful discussion with K. 
Bchnia, B. Fauque, L. Li and P. N. Ong, and correspon- 
dence with C. Ast, H. Hochst, and T. Huber. This re- 
search has been supported in part by the ICMT at UIUC 
and the NSF Grant No. DMR0605769. 



FIG. 2: (color online) (A) The Nernst signal of Ref. [D (dig- 
itized) with the surface Landau levels (n = 1, ■ ■ • 7) marked. 
The dashed lines are n = 4 and 6. (B) Fits of the inverse field 
at the peaks from Eq. (0 . The peaks identified in Ref. Q] are 
shown by circles. The straight line is obtained for g c s — 
for which BJ 1 and Bg 1 are shown by squares. The dotted 
and dashed curves are for g e s = 8 and 15, respectively. The 
quantum limit is obtained at the intersection with the vertical 
line. 



only the bulk carriers [17(. Why are the surface-state 
peaks not seen at lower fields? While we don't have 
a quantitative answer, we can qualitatively explain this 
"low- field invisibility" of surface states as follows. Firstly, 
the peaks are seen when lo c t > 1, where r is the relax- 
ation time of the carriers. Since ui c cx mT x it is smaller 
by roughly an order of magnitude at the surface than 
the bulk. The relaxation time at the surface is also ex- 
pected to be shorter than the bulk. So, the surface peaks 
can only be seen at fields higher by more than an order 
of magnitude than those for the bulk, i.e. close to the 
bulk quantum limit. Secondly, the envelope of the oscil- 
lation peaks is expected to be different in 2d and at high 
fields from the usual Lifshitz-Kosevich dependence, be- 
cause at high fields the chemical potential crosses only a 
small number of Landau levels as opposed to many in the 
low-field limit or in 3d [HI . This results in a precipitous 
decrease in the amplitude as a function of 1/-B; that is, 
the oscillations due to the surface diminish faster as the 
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